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Abstract 



We provide the first complete analysis of cosmological evolution in the Randall- 
. Sundrum model with stabilized dilaton. We give the exact expansion law for matter 

\ densities on the two branes with arbitrary equations of state. The effective four- 

' dimensional theory leads to standard cosmology at low energy. The limit of validity 

o . 

, in the high energy regime are finally discussed. 

o 

ix. 



of the low energy theory and possible deviations from the ordinary expansion law 



1 Introduction 



^ \ The past two years have witnessed the rise of several proposals for solving 

\ the hierarchy problem with the aid of large extra dimensions. In the original 

suggestion of Arkani-Hamed et al. [1], the discrepancy between the effective 
four-dimensional Planck mass and the electroweak scale originates from the 
largeness of the extra dimensions in which only gravity can propagate. How- 
ever, it is not obvious that such a large volume can appear naturally. An 
alternative scheme was proposed by Randall and Sundrum [2], and is based 
on so-called "warped compactification" . In this model, two branes are embed- 
ded into an anti-De-Sitter five-dimensional space-time (the bulk), and all the 
mass parameters of the five-dimensional action are approximately of the same 
order of magnitude. However, for moderately large values of the compactifi- 
cation radius, a strong hierarchy appears between the effective gravitational 
scale and the other mass scales in one of the branes. 

Large extra dimensions could in principle have a large impact on the cosmo- 
logical evolution of the Universe. This was first realized by Binetruy et al. [3], 
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who showed that in a two-brane model, with an empty bulk, the solution of 
the fivc-dimensional Einstein equations leads to a phenomenologically unac- 
ceptable expansion law, with the Hubble parameter proportional to the energy 
density p in the branes (rather than to ^fp as in the Priedmann equation) . This 
result changes in presence of cosmological constants [4-7]. In particular, the 
Priedmann equation is recovered at low energy, provided that the cosmological 
constants in the bulk and on the branes are related as in [2]. However, in all 
these works, a specific constraint between the energy densities of matter on 
the branes follows from the Einstein equations, and leads to negative energy 
density in our brane, at odds with phenomenology. 

Later on, Csaki et al. [8] noticed that in presence of a mechanism for the 
stabilization of the dilaton (or, more simply, of the physical distance between 
the two branes) , no constraint appears between the energy densities on the two 
branes, and the standard expansion is automatically recovered. However, in 
ref. [8] this was proved only for small energy densities constant in time. Here, 
we want to include realistic matter on the branes (with arbitrary equations of 
state) and to give a solution valid at high energy, provided that the dilaton 
is already stabilized. We are especially interested in understanding whether 
deviations from the standard expansion can occur at some energy. 

To achieve this goal, we will first give the exact solution of the five-dimensional 
Einstein equations (section I). This derivation is obtained in analogy with the 
one of [9] performed in the case of a single brane (see also [10], where some 
of our results have also been achieved). Then, in order to correctly interpret 
these results in terms of observable quantities in our brane, we will calculate in 
section 11 the gravitational action of the effective four-dimensional theory. In 
section III, we will show that at low energy the effective theory can be easily 
defined, and no deviation from standard four-dimensional cosmology can be 
observed. In particular, as in the case considered in [8], any kind of matter 
on one brane behaves exactly as dark matter on the other. Finally, in section 
IV, we will discuss the evolution of the system at high energy, and show that 
no deviation is likely to be produced when the energy density in our brane is 
smaller than (TeV)^. 



2 Exact solutions in presence of matter 

We consider the Randall-Sundrum (RS) model with two branes and the metric: 

ds^ = gAsdx^dx^ , (a,b) e {0,1,2,3,5} 

= n{y, tfdt^ - a{y, tfdx^ - bldy^ . (1) 



2 



The variable y parametrizes the extra dimension compactified on the interval 
[—1/2, 1/2]. The two branes are located at ?/ = and sX y = 1/2, and 
the Z2 symmetry y —y is imposed. We have assumed that the dilaton 
b{y, t) has already been stabilized at the value 60 • Since we are interested in 
the cosmological evolution after this has occurred, we do not enter here in 
the details of the stabilization mechanism (see however [11]). The nontrivial 
components of the Einstein tensor associated to the metric (1) are: 
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where dot denotes differentiation whith respect to t and prime with respect to 
y . We do not consider the component G55, because the corresponding Einstein 
equation accounts for the stabilization of the dilaton.Q In the RS scenario one 
introduces cosmological constants both in the bulk (A) and on the two branes 
(K)! V1/2). These energies satisfy the relation: 
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(5) 



where is the proportionality factor in the Einstein equations Gab = Tab, 
and mo a mass parameter .[^ With these choices, the system admits the solu- 
tion: 

n{y) = a{y) = e-'^"^^\y\ . (6) 



This solution gives flat space-time on the two branes. In particular, Minkowski 
space-time (with a canonical four- dimensional action) is recovered also in the 

^ One customarily assumes that a potential U (6) is generated in the 5 dimensional 
theory by some mechanism [11]. If the dilaton is very heavy, that is if near the 
minimum b we have U (b) ~ {b — b) /W' with a very high mass scale M;, , the 
55 Einstein equation gives 60 — ^ plus small corrections proportional to the energy 
densities on the two branes. If instead the dilaton is required to be fixed without 
any stabilization mechanism (that is U {h) = 0), the 55 equation forces a precise 
relation between the energy densities on the two branes, which requires the matter 
energy density on our brane to be negative [4-6]. 

^ The relation (5) is equivalent to fine-tune to zero the cosmological constant in 
ordinary Friedmann-Robertson-Walker (FRW) cosmology. 
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1/2— brane after the redefinition of the fields: 

0-^0/l]o , 1^0 = e-"^o^«/2 . (7) 

This redefinition modifies all the mass scales m of the lagrangian of the 
1/2— brane, according to 

m ^ Qnm . (8) 



This may provide a solution of the hierarchy problem if we assume that we 
live on the 1/2— brane. To see this, one first notices from the five dimensional 
gravitational action that the four dimensional effective Planck mass is here 
given by 

Mi = . (9) 

^ mo 



It thus appears natural to take all the mass scales of the system to be of the 
order of the observed Planck mass. The measured smallness of the electroweak 
scale follows then from the redefinition (8) of the masses on our brane. Since Qq 
depends exponentially on the product ttiq bo , it is sufficient to take rriQ bo ~ 70 
to account for the ratio O(10~^^) between the observed electroweak and gravi- 
tational scales. We now want to understand the behavior of the system when 
matter is included on the two branes, that is when the energy-momentum ten- 
sor of the two branes is of the form 



(^^) brane = ^ ^^"^^ (K) + Po , ^ " , - Po , K) " Po , 0) + (10) 

_^ S{y 1/2) ^.^^ / _^ ^^^^ ^ - Pi/2 , 1^1/2 - Pl/2 , Vl/2 - Pl/2 , 

Oq ^ 

where Vo,i/2 are the cosmo logical constants previously defined, while pi and 
Pi are, respectively, the density and pressure of matter on the two branes with 
equation of state Pi = Wi pi (i = 0, 1/2). In ref. [8] the system has been 
solved at first order in pi and Pi for the special case Wq = Wi/2 = —I ■ In 
the present work we provide exact solutions of the whole system for arbitrary 
equations of state on the two branes. To do this, we make use of the results 
achieved in ref. [9] in the case of a five dimensional space with one single brane. 



^ One would expect the inverse of bo to be naturally of the order of all the other 
mass scales, ^ — "i-o — • The quantity thq b^ ~ 70 needed in the RS proposal is 
however a consistent improvement with respect to the ratio between the electroweak 
and the Planck scale. Moreover, one may hope that this value can emerge naturally 
from the stabilization mechanism of the dilaton. 
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First of all, we integrate the Einstein equation (4). This is solved either for 
a = (in this case one recovers a class of static solutions including the RS 
solution (6)), or for: 

n{y,t) = X{t)a{y,t) . (11) 



This relation introduces an unknown function of time only, and considerably 
simplifies the remaining equations. Note that we have a complete freedom in 
the choice of A , since different A's correspond to different definitions of the 
time variable. By inserting eq. (11) into eq. (2), we can eliminate the time- 
derivatives, and we obtain a simple second-order differential equation for a^: 

{a'{y,t)) -4mlbla\y,t) = —^^. (12) 



This equation has a solution: 



'{t,y) = al (t) ul (y) + a?/2 (t) ^1% iv) + - ' 



2mlX {ty 



(13) 



where 
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Do 

with Co = cosh (mo &o) and Sq = sinh (mo &o) • Eq- (13) relates the value of 
a{t , y) in the whole space to the values on the two branes oq {t) = a (t , 0) 
and ai/2 it) = a{t, 1/2) . These two unknown time-dependent functions are 
determined below. Rather than the last remaining non-trivial equation asso- 
ciated to the component (3) of the Einstein tensor, we consider - in strict 
analogy to what is customarily done in conventional FRW cosmology - the 
equation associated to the identity 

= ^^^^ 



Substituting in this expression the components of the Einstein tensor with 
the corresponding components of the energy-momentum tensor, we get an 

^ To achieve this identity, the expression (11) must be used. 
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expression which is trivially satisfied in the bulk, while on the two branes it 
reduces to 



p, + 3^(A+Pi) = 0, i = 0, 1/2 
a,; 



(16) 



These two equations arc nothing but the energy-conservation law in the two 
branes and they are identical to the energy-conservation law of standard FRW 
cosmology. 

We finally have to determine the functions (t) and ai/2 (t) appearing in 
expression (13). This can be done by solving eq. (2) across the two branes. As 
shown in the work [3], one can put this last step in form of "jump conditions" 



which relate the discontinuity of n' and a' to the delta-like source (^T 
Prom the symmetry y <-> — y, we can write the "jump conditions" in the form: 
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These equations lead to the following system for oq, a 1/2: 
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As expected, the system admits no solution in absence of matter on the two 
branes, po = P1/2 — . Indeed, for this choice one recovers the static RS 
solution, which is not accounted for by the relation (11). When matter is 
instead included, the system (18) gives the solutions: 
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(19) 
(20) 
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Since A = no/ao = ^1/2/01/2, we can interpret these equations as the expan- 
sion laws of the two branes. As we will see below, eqs. (16), (19), and (20) give 
standard FRW evolution on both branes at low energy. 



3 The effective four dimensional action 



We can gain some insight on the cosmology of the RS model by integrating 
the whole action over the extra dimension y . In doing so, we make use of 
the result (13). Our goal is to get an effective four dimensional action which 
describes the evolution of the scale factors (t) , ai/2 (t) on the two branes. 

We first focus on the "purely gravitational" five dimensional action, that is we 
integrate the RS action in the absence of matter on the two branes. The latter 
will be considered eventually when we deal with the equations of motion. Our 
starting point is thus: 
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with the full (five dimensional) curvature scalar given by: 
R = 6n~2 
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Since we are interested in the evolution of the two four dimensional branes, 
we rewrite n{t , y) and a{t , y) by making use of the results of the previous 
section, eqs. (11) and (13). It is then convenient to write ^/gR and y/g in terms 
of and A: 
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With all these considerations,Q;he integral over y of the action (21) gives: 



^ The calculation can be further simplified by noticing that, from the periodicity 
imposed in the extra space, the integral of a derivative of any function of y vanishes. 
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By substituting A = no/do = ^1/2/^1/2 in the last expression^ we get 
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As we will discuss in more detail in the next section, in the low energy limit 
the equality ai/2(t) = i^oflo(^) and the related one ni/2(i) = f^o^o(^) hold. 
As a consequence, the expansion rates of the two branes are identical and the 
above action rewrites in the standard FRW form: 
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(26) 



where a = oq = f^o ^1/2 , n = no = JIq ni/2 . 

From the effective action (25) wc notice that the entire five dimensional system 
can be ( "holographically" ) expressed in terms of the physics that takes place 
on the boundaries at y = and y = 1/2 of the extra space. Notice also that 
the last term in the action (25) couples the metrics of the two walls. 

Going back to the four-dimensional action (25), and including also matter on 
the two walls, we obtain the equations of motion: 



^ In this way, we substitute A (t) with the two degrees of freedom no (t) and ni/2 (t) ■ 
The equations of motion of the effective four dimensional theory have thus to be sup- 
ported by the constraint no/ao = ^1/2/01/2 • This relation cannot be obtained from 
the action (25), since it is linked to the equation G05 = that has no counterpart 
in the four dimensional effective theory. 
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in addition to the relations which give energy conservation on the two branes, 
eqs. (16). We notice that, in the hmit po 0, pi/2 , the only solution of 
the above equations is the static RS solution ai/2 = ^odo- Moreover, one can 
verify that eqs. (27) are equivalent to the equations (19) and (20) obtained 
from the five dimensional theory. 



4 FRW evolution at low energy 

Before interpreting the four-dimensional effective theory found above, we come 
back to the static RS case. We recall that in [2] the four-dimensional metric 
g/^i, on both branes is defined as: 

g^u = n{yy^gf,^ . (28) 

The goal of this redefinition is to achieve Minkowski metric on both branes, 
in order to gain a simple physical interpretation of the system. An analogous 
procedure has to be applied also in the general case with matter on the two 
branes. 

Generally speaking, multiplying the metric by an overall function / is not 
equivalent to a change of the coordinate system. Thus, to have canonical 
normalization of the fields, the function / has to be absorbed by a redefinition 
of the fields themselves. In order to preserve the equations of motion of the 
fields, we see that we cannot choose / to depend on the coordinates t and x , 
but it can be at most a function of y . 

In analogy with what was done in the static case, we now wonder whether it 
is possible to rewrite the first four components of the five-dimensional metric 
in the form: 

9^.At,y)^f{y)9^.At) , (29) 

with g^i, of the standard FRW form diag(l, — a^, — a^, — a^). This requires the 
ratio n/a to be independent on y , that is a' /a = n'/n for every value of y . 
From the "jump conditions" (17) we see that this implies p + p = and, 
consequently, p = on the two branes. In other words, the above factorization 
is possible only if the two branes contain exclusively cosmological constants 
(in particular this is the case for the static RS solution). 
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Anyhow, it is natural to expect that condition (29) is approximately recovered 
when the matter on the two branes has a sufficiently low energy density. This 
can be understood from the results of the previous sections. From eqs. (19) 
and (20) we have: 

(1 - Ql) 

^2 77. 7W2~ 



(1 - 
12 mo 



Pl/2 



If po and pi/2 are sufficiently small, the scale factors of the two branes are 
(approximately) proportional Q by the constant factor i7o • Since n{y , t) = 
A (t) a{y , t), we have also ni/2 (t) = Qo'^o (t) ■ In particular, the ratio n/a is 
(approximately) independent on ?/ . Q 

In this low-energy limit, we can thus define the four- dimensional effective 
theory just as in the static RS model. First, we can choose the time coordinate 
so that rio and ni/2 are simultaneously time-independent. This is equivalent 
to setting X(t) = Ao/ai/2(^), where Aq is an arbitrary factor. Then, we recover 
eq.(29) with: 

/(O) = A^ /(1/2) = A^, a = X^'noao = K^a„2 . (31) 



We can now use the freedom to fix the time coordinate, and choose a particular 
value of Aq. Choosing Aq = f^o we recover, in the limit Po , P1/2 — ^ , the static 
RS solutions as presented in [2]. With this choice, the scale factor of the 
effective metrics reads a = ao = fig ^ ai/2. 

We can identify the five-dimensional quantities with those measured at low 
energy in our brane: 

• the fields must be redefined by a factor VLq. So, for instance, the observed 
density is pi/2 = Pi/2- On the other brane the canonically normalized 
density reads: po = po . 

• the total four dimensional effective action (25) acquires the form of the 
standard FRW action in terms of the scale factor a , see eq. (26). 

• The Hubble parameter of the low energy theory is given by a/a . From both 
eq. (19) and eq. (20) we get the standard evolution law: 

Notice that this relation holds exactly in the static RS case. 
^ Prom eqs. (13), (19), and (20), it is indeed possible to show that, in the low energy 
limit, the quantity n'/n — a' /a is of the same order as ai/2/ (^^o^o) — 1 • 
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while from eqs. (16) we recover 



.0, 



■p + 3-{p + p) = , 
a 

P = Po + Pi/2 , P = Po+Pi/2 ■ (33) 

Some considerations are in order. First, we would like to emphasize that at low 
energy, from the point of view of observers on both branes, the effective theory 
leads exactly to a standard four- dimensional FRW Universe. This follows from 
the fact that the standard Friedmann law is recovered, and that the energy 
densities on both branes scale with the same Hubble parameter. In particular, 
for what concerns observers on our brane, the matter on the 0-brane is re- 
garded as dark matter [8] that would completely escape any direct or indirect 
experimental detection (apart of course from its gravitational interactions). 
The gravitational effect of the matter on the 0-brane is not suppressed by 
powers of Qq , as it is the case for pi/2 . Since the only natural mass scale of 
the model is the Planck scale, po must hence be fine-tuned to small values not 
to conflict with observations (see the next section). 

Second, we remark that some care has to be paid in the interpretation of 
the physically observable quantities in the low energy effective theory. For 
instance, the alternative choice Aq = 1 [10] in eqs. (11) and (31) is not com- 
patible with the identification of pi/2 = as the observed energy density 
on our brane. This would lead to a misinterpretation of the expansion laws of 
the two branes. 

Then, in order to put quantitative limits on the validity of the low energy 
theory, we rewrite eq. (30) in terms of the observed matter densities: 

1- 



">'f^^' . (34) 

10 TeV^ 



We see that the low energy approximation is valid as long as the observed 
matter densities satisfy the bounds: 

Po < 10 M2 TeV^ , pi/2 < 10 TeV^ . (35) 



Finally, we would like to comment on Planck mass in the RS model. At low 
energy, there are two possible ways to define it, one related to the five- 
dimensional expansion, and one from the four dimensional effective action. 
These two definitions are called, respectively, local and global in ref. [12]. We 

^ We use mo ^ k~ ~ Mp and Mp Qq ~ TeV. 
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see that indeed the values of Mp obtained with these two definitions coincide 
once all the quantities in the four dimensional action are properly identified. 



5 Corrections to Standard Cosmology at high energy 



We now focus on the equations of motion when the low-energy conditions 
(35) are not fulfilled anymore. From what we said in the previous section, 
it is clear that in this regime it is not possible any longer to have a simple 
interpretation of the effective four dimensional action in terms of observable 
quantities. However, this is not important, because we make measurements 
only today, in the low-energy limit. So, it is legitimate to study the evolution 
of the system at high energy (eqs. (16), (19), and (20)), and then make contact 
with the quantities that we observe today.p^ 

We keep the previous definitions of pj, pi, Mp, and the choice A = f2o/ai/2, so 
that eqs. (20) and (16) rewrite: 



«l/2 

ai/2, 



Po + Pl/2 



12 mo 



^o) PoPl/2 
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12 mo 



Pl/2 + 3 ( Pl/2 + Pl/2 

' ai/2 ^ 







(36) 
(37) 



With our ansatz for A(t), the warp factor on our brane is constant. So, all the 
Euler-Lagrange equations on our brane are the same at high and low energy 
(i.e., they remain exactly identical to the standard equations of physics in four 
dimensions)!^. In order to close the differential system, we need an equation 



This remark should be important, for instance, when looking at cosmological 
perturbations in the early Universe. In this section, we derive only the evolution 
equations of the homogeneous background. When studying the perturbations, one 
should keep in mind that a full five-dimensional description is required at high 
energy. 

Since the freedom in choosing A(t) is equivalent to the freedom in fixing the 
time coordinate, it is obvious from general relativity principles that all physical 
results would not be affected by another choice of X{t), with the correct low-energy 
behavior A(t) Ao/ai/2(i)- It is meaningless to wonder which choice of X{t) makes 
sense physically at high energy, since contact with observations is only made at low 
energy. So, it is sufficient to give the set of equations that follows from the simplest 
choice for X{t). 
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of evolution for po- It is obtained from eqs. (16), (19), and (20): 



The differences between the evolution equations for po and pi/2 (i-e., the terms 
in the parentheses) show exphcitly that, at high energy, po is not equivalent 
to dark matter in our brane. 

Since it is assumed that mo ^ k,~ '^^^ ^ Mp and that Qq Mp ~ TeV , the above 
equations can be cast in the more transparent form : 



Ql/2 \ ^ Pl/2 (i,_PO PO \ _ PO 

ai/2j 3M2 pi/2 lOM^TeVV V lOM^TeV' 

3(po+Po) \ A _ 3(pi/2+Pl/2) \ 

2(10M2 TeV' - Po) y V 2(10 TeV^ - pi/2) J 

We now discuss the implications of these equations for the cosmological evo- 
lution in the early Universe. 

First of all, it is worth noticing that the above equations (39) encounter a 
singularity when po — lOMp TeV' . It may be possible that the presence of 
such singularity puts a limit on the theory, at least as long as the dilaton 
is assumed to be stabilized. However, as we will show later, phenomenologi- 
cal bounds from primordial nucleosynthesis indicate that this limit is hardly 
reached for pi/2 < TeV^ . 

In the regime of validity of the low-energy effective theory, po behaves as 

ordinary dark matter in our brane. So, the constraints that we usually have 
for dark matter apply to it. Although in principle we cannot say much about 
the physics on the 0-brane (in particular "non-standard" equations of state 
may be expected), we assume for simplicity that po can be decomposed into 
a constant term Pq {wq = — 1), plus matter p™ and radiation p^ components 
(with Wo = 0,1/3). 

For what concerns the constant component, the sum of the cosmological terms 
Po and is bounded by the current value of the critical density, which is of 
order 10~^'^ Mp. So, the amount of fine-tuning required here is the same as 
in usual 4-dimensional theories: 

Po + Pi/2 -P^ + ^'o Pl/2 < 10-^'^ M|. (40) 





The matter and radiation components also have to be fine-tuned to small 
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values. The best current constraint on the radiation density Pq comes from 
nucleosynthesis: since the observed abundances of light elements are only com- 
patible with an effective number of neutrinos N^ff = 3 ± 1, we see that 
is bounded by the density of one family of relativistic neutrinos. The matter 
density p!^ is obviously bounded by the value of the critical density today. 
So, in the five- dimensional theory, both p^ and p^ have to be fine-tuned to 
~ Qq p\l2 and ~ while one may naively expect po ~ P1/2 in the early 

Universe. 

Without the knowledge of the behavior of the RS model at high energy, one 

may have hoped that corrections to the standard Friedmann law could have 
solved this problem. For example, starting from po ~ P1/2 at high energy, the 
equations of motion of the system could have naturally lead to po ^ P1/2 
at temperatures of the order of the one at which primordial nucleosynthesis 
occurred. Our analysis shows that this is not the case. Indeed, let us assume 
Pi/2 ~ Po at the nucleosynthesis scale {pi ~ MeV^) and let us consider the be- 
havior of the system when it was close to the natural cut-off pi/2 ~ TeV^ . Sig- 
nificant deviations from the standard evolution are expected if at that epoch 
the energy po was almost of order TeV^ [see eq. (39)]. Going backwards 
in time, po can increase relatively to pi/2 if wq > wi/2 ■ However, assuming 
radiation domination on our brane above the nucleosynthesis scale, the above 
requirement can be met only for Wq>2, which does not seem to be a realistic 
possibility. 

A possible solution of the problem of the fine-tuning of po may arise from the 
stabilization mechanism for the dilaton, especially if it occurs at (relatively) 
low energy. Other possibilities are briefiy discussed in ref. [8]. 



6 Conclusions 

In this work we have studied the cosmological evolution of the Randall- 
Sundrum model [2] with matter on the two branes. We have first provided 
exact analytical solutions for the model, valid for arbitrary equations of state 
of the matter on the branes. By integrating the system over the extra dimen- 
sion y , we have then obtained an effective four dimensional action. 

These results can be used to investigate the physical behavior of the model. 
We have seen that at low energy the branes expand with the same rate and 
standard FRW cosmology is recovered on both of them. Prom our point of 
view, matter on the other brane is seen as dark matter. 

When one goes to higher energies, the physical interpretation of the system in 
terms of four dimensional quantities becomes less clear, since the y dependence 
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cannot be factorized away from the five dimensional metric, as occurs in the 
static case and at low energy. However, the results presented in the first part 
of this work hold also in the high energy regime. So one can still look at 
the evolution of the system at high energy, and then make contact with the 
quantities that we presently observe in the low energy theory. 

Denoting with pi/2 (respect, po) the observed energy density on our (respect, 
on the other) brane, we have found that the low effective (FRW) theory is 
valid for 

Pi/2<TeV^ , po<TeV2M2 . (41) 



As a consequence, corrections to standard cosmology cannot be found at en- 
ergies much smaller than these values. This is the main result of our paper, 
since it extends the previous analysis [8] valid only at first order in constant 
energy densities. 

The main motivation for the RS scenario is that it has only one fundamental 
scale Mp . From the definition of the physically observable quantities, this 
would suggest pi/2 ~ TeV'' and po ~ to be the most natural initial values 
for the matter densities on the two branes. While the first value appears to 
be acceptable, our analysis shows that the equation of motion of the system 
are meaningful only up to po ~ 10 TeV^ . Moreover, the phenomenological 
bound Po < Pi/2 which has to be imposed from primordial nucleosynthesis on, 
forces Po to be negligible with respect to pi/2 in the five dimensional theory. 
We have seen that the evolution of the system does not lead to this hierarchy 
at the nucleosynthesis period unless po has a non-standard equation of state. 
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